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Canonical models of single-field, slow-roll inflation do not lead to appreciable non-Gaussianity, 
unless derivative interactions of the inflaton become uncontrollably large. We propose a novel 
slow-roll scenario where scalar perturbations propagate at a subluminal speed, leading to sizeable 
equilateral non-Gaussianity, oc 1/cf , largely insensitive to the ultraviolet physics. The model 

is based on a low-energy effective theory characterized by weakly broken invariance under internal 
galileon transformations, (f> ^ (f> + , which protects the properties of perturbations from large 

quantum corrections. This provides the unique alternative to models such as DBI inflation in 
generating strongly subluminal/non-Gaussian scalar perturbations. 


Introduction: The simplest, textbook version of 
inflation consists of a single canonical scalar field—the 
inflaton—slowly rolling down a sufficiently flat potential. 
It is a common feature of these models that the magni¬ 
tude of non-Gaussianity is suppressed [ij well below the 
observable level for any foreseeable future; see, e.g., Q. 

Generically, in the low-energy inflationary effective 
field theory (EFT), there are additional higher-derivative 
corrections to the canonical action, and one may wonder 
whether these can impact non-Gaussianity in any signifi¬ 
cant way. For example, consider the following Lagrangian 


= \/^ 



V{^) + 


A4 


( 1 ) 


where V{(j)) denotes the inflaton potential and the last 
term provides a higher-order correction in the derivative 
expansion, governed by the EFT cutoff A. The virtue of 
such terms is that they do not renormalize the potential, 
while they do contribute to the non-Gaussianity Q, 

/nl ~ ^ , (2) 


eq. O can be resummed, one cannot trust values of 
non-Gaussianity greater than one. An example of such a 
resummation is provided by Dirac-Born-Infeld (DBI) in¬ 
flation Bi , where a higher-dimensional spacetime sym¬ 
metry (nonlinearly realized on (f) protects the coefficients 
of the leading derivative operators from large quantum 
corrections. For a small speed of sound of scalar pertur¬ 
bations, <C 1, DBI inflation predicts equilateral non- 
Gaussianity B , with the amplitude ^ I /c 2 BB. 

In this Letter, we propose a novel inflationary scenario, 
where the energy density of the early universe is domi¬ 
nated by the potential of a slowly-rolling scalar field, sim¬ 
ilarly to ordinary slow-roll theories. Yet, a definite set 
of higher-derivative interactions of the inflaton become 
relevant, leading to observably large non-Gaussianities. 
Nevertheless, the theory is predictive, since all the rest 
of the operators in the derivative expansion remain nat¬ 
urally small in the full quantum theory. These properties 
follow from the weakly broken invariance under inter¬ 
nal galileon transformations (9|, 

4> ^ 4> b^x^ , (3) 


where (fo denotes the background expectation value of 4> 
and /nl is the nonlinearity parameter of scalar pertur¬ 
bations [ 2 ^. It is clear from eq. m that having /nl > I 
in this model implies going beyond the (at least appar¬ 
ent) regime of validity of the low-energy effective theory, 
which requires {dcf)^ ^ A^. Therefore, unless the in¬ 
finite number of derivative operators in the ellipses of 
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that defines the Lagrangian our theory is based on. While 
the symmetry ([31) has appeared in a variety of physical 
contexts, ranging from modified gravity |9l4l2l| to scat¬ 
tering amplitudes [l^, here we use it as a guideline for 
constructing largely UV-insensitive models of the early 
universe. To our knowledge, dSD is the only alterna¬ 
tive to DBI-like symmetries for protecting the coefficients 
of strongly-coupled higher-derivative operators against 
large quantum corrections. 

Just like in DBI inflation, enhanced scalar non- 
Gaussianity is generically associated with a reduced 
speed of sound of perturbations in our model; however. 
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the enhancement is much stronger compared to the DBI 
case, the amplitude of equilateral non-Gaussianity grow¬ 
ing as /nl'*' oc 1 /cf for small c^. 

The model: The theory we wish to study below is 
defined as a combination of an inflationary potential and 
the four Lagrangian terms Q , which, barring the factors 
of v^— 5 , can be written as follows 


L2 = K\G2{X), (4) 

L, = AiG3iXm, (5) 

L4 = M^^G4iX)R + 2A\G'^{X) ([$]2 - [$2]) , (6) 

- \aIG'^{X) ([$]3 - 3[<i>][<l>2] + 2[$3]) . (7) 


Here, $ is a matrix, consisting of second derivatives of 
the inflaton, and the brackets [...] 

denote the trace operator. Moreover, Ga are arbitrary 
dimensionless functions of the dimensionless variablqj 
X = A 2 ■ For simplicity, we will assume 

that these functions can be Taylor-expanded around zero, 
Ga{X) = Furthermore, the two scales in 

the theory are related to each other as A 2 = MpiA|, so 
that the smaller of these, A3, can be regarded as the gen¬ 
uine cutoff of the underlying low-energy EFT (we assume 
Mp\ is the parametrically highest scale in the problem). 

In the limit Mpi —>■ 00 , the Lagrangian terms in Q- 
© reduce to the galileons of Ref. Q , which are exactly 
invariant (up to a total derivative) under ([3]). In this 
limit, there is a non-renormalization theorem, accord¬ 
ing to which the galileon operators are not corrected, 
at least perturbatively, by quantum loops 1^. For a fi¬ 
nite Planck mass, the operators (01)-© break the galileon 
symmetry, but only weakly Q. The defining property of 
the theories with weakly broken invariance under dS]) is 
a generalization of the non-renormalization theorem of 
Ref. 1^, which renders the quantum corrections to the 
coefficients suppressed by positive (integer) powers 
of the tiny ratio Aa/Mpi Q. In addition, these theories 
belong to the so-called Horndeski class of scalar-tensor 
models [l^, characterized by second-order equations of 
motion both for the scalar and the metric fl6| . 

It has been shown in Ref. Q that the properties of 
the theories with weakly broken galileon symmetry imply 
the possibility of a moderately coupled, yet predictive, 
regime characterized by 


(bn 




( 8 ) 


for a homogeneous (/(-profile on a FRW background with 
the Hubble rate H. From now on, X will be under- 


^ We use the ‘mostly plus’ signature for the metric. 


stood as evaluated on the background solution. Despite 
the moderate coupling, quantum corrections are under 
control when the scalar background profile satisfies ©— 
even in the case that these inequalities are saturated— 
and the predictions of the classical theory can be trusted. 

The slow-roll backgrounds: As noted above, we 
will be interested in the potentially-dominated models 
of inflation, characterized by weakly broken invariance 
under the galileon transformations. These are governed 
by the following Lagrangian 


^ 


A/f 2 1 

—- - ((9(/>)^ - H((/)) + y^ Z/t -h ... 

i=2 


, (9) 


where, since we have extracted the canonical scalar and 
graviton kinetic terms, G 2 is assumed to start at least 
quadratic in X, while G 3 can have a linear piece. From 
now on we will set G 4 = G 5 = 0 for the sake of simplic¬ 
ity; generalization to the case of nonzero G 4 and G 5 is 
straightforward and will be commented on where appro¬ 
priate. The ellipses in © denote an infinite number of 
other operators, present in the low-energy effective the¬ 
ory. We will assume that the potential V{(j)) satisfies 
the ordinary slow-roll conditions, ey <C 1 and \riv\ 1 , 
where the (potential-based) slow-roll parameters ey and 
rjv are defined as 


ey = 




r]v = Alpy 


XI 

V 


( 10 ) 


The previous analysis of quantum loops, leading to 
the non-renormalization theorem summarized above, has 
concentrated on the case with a vanishing potential [^. 
It is straightforward to show that the same results remain 
intact also in the presence of a nonzero, but sufficiently 
flat V((/>), satisfying (ITOl) . 

For the flat FRW ansatz, ds^ = —dt^ + a?{t)d3p‘^ the 
two Friedmann equations that follow from © read 


= V - AiX 


Tl 

f 2 




2A/|iiZ = -AjXF{X, Z) -H 2MpWG30o 


( 11 ) 

( 12 ) 


where H = d/a and the function F{X, Z) = 1 + 2 G 2 — 
dZG'^ has been introduced for later convenience. More¬ 
over, in the slow-roll regime the homogeneous equation 
of motion of (p reduces to 

3HPoFiX,Z)c:,-V\Po) ■ (13) 

We are interested in a regime where higher-derivative op¬ 
erators in © become important, while the quantum cor¬ 
rections are still under control. To this end, we assume 
Z ^ 1 , which also fixes the magnitude of the parameter 
X. Indeed, from the definition of X and Z, eq. ©, it 
follows that = A|Z/(MpiiZ); making use of eq. m, 
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one immediately obtains X ~ ^/e, where e = 

Note the order-unity slowly varying function of time 
F{X^ Z) in eqs. (fTT]) -(fT 3 )) . which is strictly one in canon¬ 
ical slow-roll inflation. Apart from this minor modifi¬ 
cation, all the equations that describe the background 
solution are similar to those of ordinary slow-roll models 
(up to corrections of higher order in ty and rjv)- In par¬ 
ticular, unlike e.g. the DBI case, the usual flatness con¬ 
ditions ey <C 1 and \r]v \ 1 need to hold for sustaining 

the quasi-de Sitter phase in our model. It is precisely for 
this reason that we refer to it as “slow-roll” @ At the level 
of perturbations, our scenario is of course very different 
from the canonical slow-roll inflation; for example, unlike 
the latter, the scalar perturbations become strongly cou¬ 
pled at an energy scale parametrically smaller than Mpi, 
something we discuss in greater detail below. 


It follows from the Friedmann equations that the con¬ 
tributions from the derivative operators in dO]) to the in¬ 
flationary energy density and pressure are proportional to 
X ^ Y^. One may wonder therefore, whether loop cor¬ 
rections can outweigh these contributions for small values 
of e. For Z ^ 1 , the leading quantum corrections to the 
background stress tensor scale as ~ A| Q. This should 
be much smaller than which implies a lower bound 

on the slow-roll parameter, e ^ [H/My\f'. This is the 
same bound on e as the one that arises from requiring 
quantum fluctuations of the inflaton to be small 22, 2^. 

Non-Gaussianity: Inflationary theories can be 
conveniently studied in a model-independent way using 
the EFT framework 1 ^, 2 ^. To this end, we decompose 
the metric in the ADM variables, 


ds^ = -N^df + 7y {N^dt + dx^) {N^dt + dX) , ( 14 ) 


and work in the unitary gauge, where the constant¬ 
time hypersurfaces are chosen to coincide with those 
of uniform (f). The perturbed quantities are defined as 
6N = N — 1 , and SK = K — 3 H, where K denotes the 
trace of the extrinsic curvature of equal-time hypersur¬ 
faces. The action in eq. dH]) (with G4 = G5 = 0 ) can 
be expanded to the cubic order in perturbations in the 
following way 211 


S= d'^Xy/^ 


- M 2 i( 3 ij 2 + H) 


Ar2 




6N^ -f M^SN^ - Ml5K5N + M^SK6N^ 


( 15 ) 


^ Inflationary models based on particular subsets of the Lagrangian 
terms in @-0 have been studied in Refs. [T^ . 11 ^ . However, 
these references have focused on kinetically-driven inflation, cor¬ 
responding to A 2 ~ Mp^H^ and X Z ~ 1 in our notation. 


where M2, M3, Mi and M2 are functions of time of 
canonical dimension one. In terms of the functions G2 
and Gs introduced in eqs. o and m, these are 


= -2A^A 
= -2A^A 


3ZG3 + 6ZXG3 - 2XG2 - YG's 
3 XGo + -A^G'," 

43 2 


- Z (4G3 -f IIAG3 -h 2X^G3) + YG's 


M! = -2^ZG' , 


= -2^^ (2G^ 


AG'' 


( 16 ) 


where Y = (/(g/Ag. One can see from (ITS)) , that approx¬ 
imate invariance under galileon transformations imposes 
the following (radiatively stable) hierarchy among the 
various EFT coefficients: 

^ ~ M^,H, Ml ~ ^ . ( 17 ) 

This is in stark contrast to what happens e.g. in solely 
shift-symmetric theories, where the coefficients that stem 
from higher-derivative operators such as Mf and M| i 
much stronger suppressed. The latter hierarchy moti¬ 
vates to define the dimensionless, order-one coefficients 


0:1,3 = ~ 


M |,3 
2A/2j77 ’ 


02,4 = — 


MI2H 

2M2ji7 ’ 


( 18 ) 


convenient for describing the parameter space of the the¬ 
ories at hand. 

At sufficiently high energies (encompassing the Hubble 
scale), the dynamics of scalar perturbations is fully dom¬ 
inated by the dynamics of the adiabatic mode tt, defined 
through (j){x,t) = (j)o {t + ^{xjt)) [l^, In the decou¬ 
pling limit corresponding to this regime, the scalar part 
of the action (USD reads 


(1 -I- ai) TT^ - c; 


= / d'^x X 


-I- (a2 - Oi) IT -—5- 2 (ai + 0:3)71'^ 

02—04 .2 d^TT 02 {d'K^'d'^TT 
+ 2 -77- TT ' 




H 


H 


( 19 ) 


where the speed of sound is 


2 1 -f 02 

Cs = T—- • 

1 - 1-01 


( 20 ) 


It follows from eq. ( 1201 ) that if, for whatever reason, the 
parameter 02 happens to be close to —1, one can have 
strongly subluminal scalar perturbations. Most impor- 
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tantly, the approximate galileon invariance guarantees 
that such an ‘accidental’ arrangement of the parameters 
is respected by loop corrections. This is qualitatively dif¬ 
ferent from how the small arises in models such as DBI 
inflation, as we discuss in detail belovU. 

It is worth stressing at this point that the operators in 
the last line of (US can be rewritten in terms of those 
in the second line via a perturbative field redefinition 
0. This simply amounts to using the linear equation 
of motion in dm). After straightforward manipulations, 
one flndtll 





.[dnf 


■ 727 r 


( 21 ) 


The two operators in (|2T]) are precisely those appearing 
in the decoupling limit of DBI theories and the bispec¬ 
trum they produce is close to the equilateral shape 0. 
The genuine difference arises once the magnitude of non- 


Gaussianities is concerned: instead of the 


1/c? 

behaviour characteristic of DBI inflation [ 7 [, in theories 
with weakly broken galileon symmetry non-Gaussianity 
scales as ~ l/c^ in the small-Cg limit. The latter 

scaling is due to the last operator in (IT^ . whose precise 
contribution to the three-point function of the curvature 
perturbation C, reads 18 , lin 


B(^{ki,k2,kz) 


1 02 1 k\{k\ — k 2 — k^) 

16 1 + ai kt{kik2kzY 


X 1-H 


^iyjkikj ^kik2k^ 




-b3- 


k! 


+ 2 perms . 


( 22 ) 


Here, A denotes the normalization of the 7r-kinetic term 
in the decoupling limit, A = (—MpjiJ)(l -|-ai). The am¬ 
plitude of non-Gaussianity can be directly read off from 
eq. (1^ 

.equii _ 5 k^B(;{k^,k^,k^) _ 65 0:2 1 

/nl - ^ 2 ^ ~ 162 1 +aicy ^ ’ 


^ In the limit X ^ 0 , using eq. (| 16 l> with G4 = G5 = 0 in eq. 
one finds ai = 802 • This gives a negative kinetic term to tt for 
02 — —I- The parameter X need not be very small, however; it 
is of order X \fk a few x 0.1 in e.g. slow-roll models with 
monomial potentials. Moreover, the relation oi = 80:2 no longer 
holds for G4 7^ 0 or G5 7^ 0. 

Explicitly, the coefficients 71 and 72 read 

71 = (c?-l) (1 + ^) +(l + c2)«i , 

72 = 2 ^ ^2 -I- —^ -I- ^(“1 — 204) + 2 ai — 203 . 

^ We follow the standard definition of the comoving curvature 
perturbation, gij = a^e^^Sij. The three-point function is de- 

fined as {({ki)({k 2 )({k 3 )) = (j 2 iki^ B^{ki,k 2 ,k 3 ), 

and kt = ki k2 ks . 


where = fc^P^(A:) = H^/{ 4 Acl) is the dimensionless 
power spectrum, evaluated at a fiducial momentum scale 
/c*. A significantly reduced speed of sound, 02 — —1 
(see eq. (120)1 1. implies a negative /nl • However, due to 
the strong dependence of /nl on c^, even slightly sub¬ 
luminal perturbations can produce a sizeable amount of 
non-Gaussianity; for example, a 10 % tuning of the 0:2 
parameter can give rise to /nl'^'* — ~70 for ai = 1 and 
02 = — 0 . 9 . We stress that such a tuning is not ‘un¬ 
natural’ as a result of the non-renormalization theorem 
outlined above. 

If the theory (fT^ is to be predictive, it is crucial that tt 
is weakly coupled at energies of the order of the inflation¬ 
ary Hubble rate. A* ^ H, where A,^ is the energy scale 
at which perturbative unitarity is violated in the 2 —>■ 2 
scattering of tt. In the <C 1 limit. A* is set by the 
last interaction term in the action ( 1191 ) . and is estimated 
as A* ^ AaCs^^^. For oi ~ 1 and 02 — — 1 , using the 
experimental value ~ 57r^ x 10“®, one finds 

A’ ~ 1= mf. ( 24 ) 

|/nl I 

Even for the largest /^ 7 ’^ compatible with the current 
observational bounds [^, the strong coupling scale is 
fairly above H, but well below the symmetry breaking 
scale 0, A 3 ^ 0 ( 5 )|/X> 2 - 

We close this section with a remark concerning the 
regime of validity of the decoupling limit. For the small 
values of the speed of sound we are interested in, one 
should be careful with mixing terms that involve spatial 
derivatives. For example, consider the 6 NSK operator 
in eq. disi). The most important mixing of scalar modes 
with gravity that arises from this operator is suppressed 
by a factor of e/c^ compared to the 7 r-kinetic term at 
horizon crossing. Therefore, the validity of the decou¬ 
pling limit analysis requires that cf > e hold, which puts 
an upper bound (/nl ^ 1/e^) on the amplitude of non- 
Gaussianity attainable within the decoupling limit. 

Discussion: In ‘ordinary’ inflationary theories, the 
statistical properties of perturbations are mostly deter¬ 
mined by operators with the least number of derivatives, 
i.e. ( 9 (/))^” or, in the EFT language, by operators of the 
form (SN)'^. The prototype example of this is provided 
by DBI inflation, where the background can be consis¬ 
tently strongly coupled, with an infinite number of oper¬ 
ators of the above type becoming relevant for its dynam¬ 
ics. In the EFT of eq. 013 , this corresponds to the large 
coefficient or, equivalently, to oi 1 in the nota¬ 
tion of eq. dUD, and implies a parametrically suppressed 
speed of sound, cl ~ I/oi, see eq. (l20l) (higher-derivative 
operators are effectively negligible in DBI inflation, so 
one can set 02 — 0 ). Despite the strong coupling, the 
symmetries of DBI theories protect the structure of the 
Lagrangian from large quantum corrections. This mech- 
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anism of obtaining strongly sub-luminal scalar perturba¬ 
tions has been extensively studied in the literature (see, 
e.g., and provides an attractive way of gen¬ 

erating large equilateral non-Gaussianity in single-field 
inflation. 

In this work we have proposed an alternative scenario 
that allows for strongly subluminal scalar perturbations 
within a well-defined low-energy EFT. Our model re¬ 
lies on moderate coupling, i.e. ~ 1, as a result of 
which more than one operator in eq. m become large 
enough to affect the scalar perturbations significantly. In 
the simplest realization considered above, the scalar two- 
and three-point functions are determined by the opera¬ 
tors and 6N6K. Instead of the ai ^ 1, 012 “C I 
case characteristic of DBI inflation, our mechanism relies 
on an adjustment of the order-unity ai and a 2 , which 
results in a somewhat reduced speed of sound in (l20l) . 
The theories we have studied in this work not only al¬ 
low for such an adjustment, but also provide a way to 
protect it against loop corrections. The central reason 
behind the robustness of the results obtained in the clas¬ 
sical theory is the weakly broken galileon symmetry Q, 
which inherits the remarkable quantum properties of the 
galileon operators Q, exactly invariant under ([3]). 
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